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Abstract 


This  paper  consists  of  two  notes  on  the  equations  governing  the 
propagation  of  electromagnetic  wares  in  the  ionosphere.  In  Part  I  the  ten- 
sor properties  of  the  ionosphere  are  derived  by  a  method  due  to  van  der  Wyck 
(tlj  »  PM  ).  While  we  introduce  only  minor  improvements  on  his  method,  it  v;es 
felt  worthwhile  to  make  these  results  generally  more  available* 

Pert  II  shows  the  equivalence  of  various  forms  of  the  propagation 
equations.  A  general  scheme  is  indicated  for  finding  simplifying  transforma- 
tions of  the  equations  in  their  first-order  from.   In  particular,  we  bttain 
that  transformation  which  yields  Rydbeck's  coupled  equations,  and  give  an 
interpretation. 
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PAEI  I 
nSE  lENSOR  FROPERIIES  OF  THE  lONOSPHEEE. 

1.   Introduction 

It  is  well  knovm  that  the  linearized  macroscopic  field  eq^iations 
in  the  ionosphere  involve  tensor  properties  of  the  medium.  The  most  accurate 
derivation  of  these  equations  shoiild  start  with  the  microscopic  field  equations, 
and  apply  the  methods  of  statistical  mechanics  to  yield  the  tensor  field  equa- 
tions for  the  macroscopic  fields  (see[l]).  Some  such  attempts  have  "been  made 
by  Hartree  [2]  and  Darwin  [3]  for  a  case  in  which  the  ionosphere  may  he  describ- 
ed by  scalar  properties.  However,  the  general  derivation  in  which  tensor  proper- 
ties arise  seems  quite  difficult  and  remains  to  be  carried  out.  In  place  of  this 
we  are,  at  present,  forced  to  accept  derivations  in  which  the  averaging  is  less 
detailed.  Many  such  derivations  have  been  given  but  they  are  not  as  rigorous  as 
might  be  expected,  and  contain  numerous  assucrptions  which  are  not  explicitly  stat- 
ed. In  the  present  note  we  present  a  modification  of  a  derivation  due  to  van  der 
Wyck  [4J  in  which  the  averaging  is  explicit.   The  main  assumptions  are  stated  and 
it  is  shown  that  the  resultant  tensor  properties  of  the  ionosphere  (when  some  of 
van  der  Wyck's  simplifying  assumptions  are  eliminated)  are  just  those  obtained  by 
the  less  rigorous  methods,  dis  may  serve  as  a  justification  for  the  use  of  the 
usual  derivations.  Por  comparison,  one  of  the  standard  derivations,  essentially 
that  of  Baner jea  [5]  ,  is  included. 

2,   Tormulation 

The  general  macroscopic  field  equations  are 

^7K  E  -  -  ^^ 


(1) 


o'C 


together  with  the  relations 

(2)  ?  =  c^  +  ? 

T  =  0-?+  t . 
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Here  jul  ,c  and  cr  have  their  usual  meanings.  P,  the  polarization,  is  a  measure 
of  the  electrical  "  stress"  of  the  medium, which  is  proportional  to  the  average 
displ£icement  of  the  electrons  from  their  equilibrium  positions,  Z^   the  total 
current  density,  consists  of  the  usual  conduction  current  density,  cTE,  and  a  con- 
vection ctirrent  density,  I,  which  is  the  average  flow  of  free  charge  density. 
The  ionosphere  is  represented  "by  free  space,  i.e.,  e  = /t-  =  1  and 
o^  =  0,  in  which  there  is  a  given  density  distribution,  N(x,y,z),  of  free 
electrons.  Heavier  particles (ions)  are  assumed  to  be  present,  and  we  assvune 
their  effects  may  be  taken  into  account  by  means  of  the  considerations  given 
in  Section  h,     A  static  magnetic  field,  H_,,  due  to  the  earth,  is  also  present. 

The  problem   is  to  derive  expressions  for  the  macroscopic  quanti- 
ties P  and  I  in  terms  of  the  field  conrponents  and  ionospheric  properties.  The 
equations  of  motion  for  charged  particles  in  the  presence  of  electromagnetic 
fields  mast  play  a  major  part  in  this  derivation;  these  equations  are  introduc- 
ed below. 

?,   Standard  Derivafcion 

In  the  standard  derivation  it  is  assumed  that  1=0  and  that 

(3)  ?  =  N  e  ?, 

where  p  is  the  average  displaxiement  of  an  electron  from  its  equilibrium  position 
and  e  is  the  charge  on  an  electron.  To  coOpute  p  one  considers  the  equation  of 
motion  for  an  '•  average"  electron: 

C)     ?*"?  -  S[?-%]  =  l^  . 

where  m  is  the  mass  of  an  electron,  and  where  we  have  used  the  linearizing  ap- 
proximation Jh]  <<  l^rjl.  The  electric  force  term,  —  S,  is  taken  here  to  be  the 

valioe  derived  by  Darwin  [3J  and  generally  used  in  ionospheric  problems,  although 
this  has  been  questioned  by  some  authors.  The  so-called  friction  term,  i'  p  , 
arises  from  the  fact  that  iS)   is  the  equation  of  motion  for  an  average  electron. 
It  is  assumed  that  collisions  of  the  electrons  with  the  heavier  ions  occur  with 
average  frequency  "^ .     The  effect  of  these  collisions  is  assumed  to  be  a  viscous 
damping,  as  is  indicated  by  the  equation.  However,  the  presence  of  this  term  is 
one  of  the  main  points  which  require  a  more  detailed  averaging  process.  In 
standard  derivations,  such  as  that  of  the  present  Section,  the  friction  term  is 
assumed  and  hence  eo^uation  (^)  is  used,  while  more  detailed  considerations  such 


-> 


as  those  of  Darwin  inTolve  the  complicated  averaging  process  of  statistical 
meclianics. 

If  we  assume  p,  E,  and  H  to  have  time  factors       of  the  form  e  "'  , 


eqiiation   (4)  "becomes 


(5) 


(l-lZ)?+i   [?.?i     =    .^i'    . 


Ne 


where  we  have  introduced  the  notation 


mw 


(6) 


^         TTBIl  • 


1  -  i  Z     . 


In  matrix  form,  (5)  is 


(7) 


P2 

-X 

He 

\    -=^2    S'l     "  / 

1^3^ 

\^/ 

Inverting  the  coefficient  matrix  yields  the  average  displscepent 


(8) 


Ne6(6^-  Y^) 


^1  "  ^^    ^1^2  "*"   ^^^3 


\  ^3^1  *  i«i^2  =^3^2  -  -^h 


7^7 J  -  1672  1 


72^3 


+  i6y. 


2   .2 

73-6     y 


f'M 


\  E, 


-  Ne 


1  -^  . 


Here  we  have  introduced  the  3x3  matrix  T, 

2     2 
The  ahove  procedure  is  invalid  when  6  =  Y  =  Ij  this  is_ji  resonance 

which  occurs  when  the  wave  freq[uenc7  tu  equals  the  gyrof re  quency  I  ^TI  |  and  when 

there  is  no  damping. 


If  the  opposite  time  factor,  e~   ,  )isd  been  used,  the  only  effect  would  "be 
to  -change  the  sign  of  i  throughout  the  remainder  of  Pert  1.  The  resulting 
ejcpresnions  are  given  in  Part  II  where  this  tine-dependence  is  assumed. 
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S\ibstltutine  (8)  in   (3)  we  have 

P     =     T  ?    , 
and  Maxwell's  eqimtions    (l)  tecone 

7   X  E    =  -  iu)  H 

(9)  ^ 

7   X  H    =     io)  K  E 

where  we  have  introduced  the  ionospheric  dielectric  tensor 

(10)  1  =  T  +  T  . 

Here  I  is  the  unit  matrix. 

Tor  purposes  of  coniparison  with  the  result  of  the  next  Section  we  may 


take  y,  =  yj,  =  0,  y«  =  y  (i.e.,  we  take  £,  in  the  vertical  direction),  in 
which  case  T  "becomes 


(11)       T  = 


(62-y2) 


-8 

iy 

0 

-iy 

-6 

0 

0 

0 

6 

/ 


4.   Hiporo-'as  Derivation_  (see-  van  der  Wyck:.[4j.  p.  IQ-IQ) 
Let  P  =  0  and  take 

(12)  f  =  He  ^  , 

where  v  is  the  average  velocity  of  electrons  due  to  external  electroma^-netic  field 

(i.e.,  a  perturTsation  on  the  average  thermal  velocity  V  ),  To  compute  v  we  take 

__  e 

the  coordinates  such  that  H_,  =  (0,0,H_)  and  consider  the  equation  of  motion  for 
an  individxie.l  electron  "between  collisions: 

(«)  ^-  SC"-^  -  i^  • 

Here  v  is  the  velocity  of  the  electron  due  to  the  external  fields.  We  introduce 
the  Quantities 


(14) 


v+  =  v^  i  iv^ 


\    -    \    t    ^2 


-5- 
into  (13)  aad  obtain^ componentwise, 


(15) 


iwy  T^ 

-     m\ 

• 

r      - 

iup-  T_ 

• 

^3 

=     1^3 

Here,  as  in  (6),  vre  have  used  y  =  ,  iVom  the  esponential  time  dependence 

—    mu) 

of  E,  as  ess-umed  in  Section  3»  w®  °ust  hare 

(16)  E.     =    B^  8in(«jt  +  f)_^) 

where  ^  are  phase  constants.  Thus  the  first  of  equations  (15) »  which  is  a 
simple  first-order  linear  inhomogeneous  equation,  has  the  general  solution 

(17)  v^    =  e-^y^*  I  J^+/e^^"'*  si^(«'t  +  ^+)  "it  +  co^st  j  . 

To  evaluate  the  constant  it  is  assvuaed  that  |V  [  »  |v| ,  i.p.,  that  the  theraial 
velocity  is  much  larger  than  the  perturbation  velocity;  hence  the  velocity  v 
after  a  collision  is  taken  to  be  zero.  If  the  last  collision  occurred  'X   seconds 
before  the  time  t  the  solution  (l?)  subject  to  the  above  assuiiiption  becomes 

cos  («,frf^^ )  [  (e-^^  ^  cos  u,  t  -1)  +iy  (  e-^^Ui  ^  ^^^  ^^'\ 

(18)  P  n"^ 

+  8in(«,t+^^)|  (e-^^ru-r  ^^^  ^^)  _iy  (^-lyu^r  ^^^  ^^  ^^ 

This  solution  vanishes  when  X    =  0,   i.e.,  where  a  collision  has  occurred  at  time 
t.  "So^-r   it  is  assumed  that  the  probability  that  the  electron  had  a  collision  in 

interval  dT  around  t   is  given  by 

V       Y 
(19)     f-  exp(-  -fr)  ^X     =  V  esp(-yr  )  dT 

where  V  =  average  thermal  velocity,  J^  =  mean  free  path  of  the  electrons^  and 

T 
V    -  0  is  the  average  collisional  frequency.  The  average  velocity  Y  is  then 

"    r 


r 
*      "  (i-y^) 
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(20) 


ir^(t)  =/°°  T^(t.r) 


V  e 


-rr 


dr 


and  from  (18)  we  have 


^21)  \=  i£i|-«o«(«'*^^+)  i^(afi7)^  ]  "  '^"^^'*^-^^ 


(Zfly) 
l+(Zfiy)2j  ^ 


Similarly  we  obtain 
(22)   V  =  -2-S  - 


mu)  - 


-cos 


(uiW-^  ) r  I  +  sin(wt+^  )| 

-Li+(z-iy)^J         "L 


(2^iy) 
l+(Z-iy)2j  ^ 


Uow,  recalling  that  ^  E,  =  iw  E,,  we  derive  from  (l4),(2l),  and  (22) 

ot  J      J 


(23) 


"^  (fiV) 


-  16 


Trom  the  last  of  equations  (15)  we  readily  obtain 

Equations  (23)  and  (2^)  s-uistituted  in  (l2)  give  the  current  density: 


T      =     -L    3i_E 
"3         ffiuj      6    ^3    • 


^\ 


(25) 


iujX 


(sV) 


-6 

0 


0 

0 


=     ioj 


S    ^         e 


Maxwell's  eq-uations    (l)  again  reduce  to    (9)   if  we  introduce 


(26) 


E    =     I  +  S     . 


From  equations  (lO), (ll), (25)  and  (26)  we  see  that  the  dielectric  tensors  agree  for 
the  present  choice  of  coordinates.  By  rotating  these  coordinates  to  any  arhitrary 
position  it  is  easily  seen  that  they  must  agree  in  general. 

The  derivation  of  the  present  Section  is  that  of  van  der  Wyck:£4j  up  to 
equations  (21)  and  (22).  However,  in  going  to  (23)  he  introduces  some  asstanptions 

which  are  apparently  unnecessary  and  make  his  expressions  different  from  those 
usnally  foxmd  in  the  literature. 
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PART  II 
EQtnVAlENT  rORMS  0?  MAXWELL'S  EogATICBTS  FOR  SHE  lONOSPHEEES 

1.  Introduction 

Maxwell's  equations  in  the  ionosphere  have  "been  considered  "by  various 
authors,  and  many  of  these  investigations  use  different  forms  of  the  equations  for 
similar  configurations  of  incidence  and  earth's  magnetic  field.  The  reasons  for 
using  different  "but  eqiiivalent  forms  of  the  equations  are  varied  but  generally  in- 
volve some  atteinpt  to  reduce  the  cotgjling  "between  ordinary  and  extraordinary  waves. 
In  the  present  note  we  show  a  general  procedure  for  obtaining  different  forms  of  the 
field  equations.  Prom  the  method  outlined  it  is  easy  to  see  the  connection  between 
some  well-known  forms.  In  particular  the  equations  for  the  n  -  waves  of  Eydbeck[lJ 
are  deduced  and  the  quantities  tt.  and  tt  are  shown  to  come  from  a  rotation  of  the 
field  quantities  throxigh  a  complex  angle,  A  similEir  treatment  for  a  different 
ionospheric  configuration  is  considered,  and  the  results  seem  to  be  unnecessarily 
complex.  This  is  an  indication  that  the  analysis  of  Itydbeck  is  not  readily  ap- 
plicable to  configurations  other  than  the  one  he  considered, 

2,  G-enersl  rormolation 

For  electromagnetic  fields  in  the  ionosphere,  with  time  dependence 
of  the  form  e~   ,  Maxwell's  eqimtions  are 

U  X  f  =  i  i!if 
c 

(2.1)  .       ,  _      . 

V   ^   f   =     -i|z(z)E  . 

The  ionosphere  is  assumed  to  be  horizontally  stratified,  and 

1-A(62_y2)     ACy^y^+iSy^)     ACy^y^-iSy^)  \ 


(2.2)  KT)     =  (k^j(z))  =  I  ACy^y^-iey^)     1-A(6%'|)  ACy^y^+iSy^) 

ACy^y^+ifiy^)    A(y2y3-i6y^)  l-A(6^-y|)   / 

Equations  (2,l)  and  (2,2)  are  derived  in  Part  I.  The  notation  in  equation  (2.2)  is 
related  to  the  more  standard  symbols  as  follows: 


(2.3) 


-9- 


fi    =    1  -  i  Z(z) 
Z(z)     =    ^^ 

'         •—  111 


A     =     X(2) 


(l-iZ)[(l-iZ)2»T^ 


ui 


^'C.)   =    (71.y2.y3)   ^  =^  ^(^) 


ZU)     =     ^N(z) 


sxu 


We  seek  solutions  of  equations    (2.I)  in  the  form  (see  [2]) 
^(:^.y.z)     =     ^(z)  esp(  i  I  [px  +  qy]) 
^(x.y.z)     =     I^z)  exp(  i  ^  1^2+  qy])        . 
When  equations    (2.4)  are  substituted  in   (2.1)  we  have 


(2.4) 


r 


(2.5) 

and 
(2.6) 


v^^  =  fc  K^^^'-i^i^^^ -^  Vi(-)  *  k-oV^a 


33 


1^     13(3)     =      |i)72(z)  -  qVi(z)] 
f^    =     i^A(z)W(z)      . 


Here  we  have  introduced  the  matriz  notation 


IC33         iC33 


-q 


33 


(2.7)     W(2)  = 


l2(.) 


1^(2) 


J   A(z)  = 


,^^2,^11^33-^13^31]  _^    ^  ^2 
^33  J         33         ^33 


^33 


pqt 


'^12^33-^3'^32 


^33 


^33        ^33 


2_^2^2li2i22 


33 


M  ...       ^        \    ^33      ^33     ^33        ^33 

Maxwell's  equations  for  an  arbitrarily  incident  plane  wave  in  the  ionosphere  have 
thus  heen  reduced  to  the  first-order  systoEi. 
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Howerer,  it  is  freqiiently  advantageous  to  introduce  new  dependent 
varialsles  which  are  linear  comoinations  of  the  7^^,  V  ,  I^,  and  I  .  Thus  if 
P(z)  is  some  non-singular  4^4  matrix  we  may  define  the  new  variahle  "by 

(2.8)       UCz)  =  P"^(2)W(z)    or    W(z)  =  P(z)U(z) 

The  system  of  equations  satisfied  "by  these  new  variahles  is,  from  equations  (2.6) 
and  (2.8), 


(2.9) 


dTJ(: 


dz 


■1.  [i^p-i«.r^g 


1U)  . 


The  motivation  in  choosing  a  transformation  matrix  P(z)  is  usually  to  put  equation 
(2.9)  into  a  form  which  can  he  readily  solved  (or  solved  approximately;  see  l2J»l3j)< 
However,  in  the  present  note  we  shall  use  this  method  to  obtain  and  compare  some 
equivalent  forms  of  Maxwell's  equations  for  the  ionosphere, 

3.  Uormal  Incidence  (Oblique  magnetic  field) 

For  a  plane  wave  normally  incident,  p  =  q  =  0  and  the  xy  -  plane  may  "be 
rotated  ahout  the  z  -  axis  so  that  the  earth's  magnetic  field  £,  shall  lie  in  the 
yz  -  plane,  say.  Therefore^in  this  case  we  may  take  y^  =  0  in  eq.  (2.3)  without 
loss  of  generality.  The  coefficient  matrix  of  the  system  (2.6)  hecomes 


A(z)  = 


where 

^a  =  1  -  X  8  ii-^ 


(2.11)  <P  =  [(82-  y2)  -  xe]  -i^^i^  ,      D  =   r(6-X)(82-T2)  -xy^j   . 

:  ^  ixy3li^ 

The  term  V  serves  as  a  measure  of  the  coupling  between  "  waves  •' determined  "by  the 
orthogonal  pairs  T.,  I  and  I-,  V  of  the  field  components.  Thus  if  V  were  zero 


-li- 


the system  (2,6)  would  reduce  to  two  independent  systems  of  two  equations  each. 
This  is  readily  seen  if  we  eliminate  I.  and  I.  from  the  system  (2,6),  using  the 
coefficient  matrix  A(z)  given  hy  (2.10),  We  then  obtain  the  equations 


dV(z)      2  2 

\-    +  ^    a(z)T^(z)  =  -^y{z)l^{z) 


(2.12) 


dz 


d^„(z)     2  2 

-%-  +  ^  P(z)  yz)  =    J^Y(z)V^(z)  , 


dz' 


which  uncouple  if  V  =  0,  Ihe  vertical  field  components  are,  from  equations  (2.5), 


13(2)  =  0 

-3Sy3 

3'   "  [(6-X)(62-y2).6y|l 


ri8V^(z) 


^7^^^U)]       . 


These  equations  have  heen  obtained  and  studied  by  Ginsberg  \k\   and  otiiers.  The 
coupling  coefficient  enters  into  equations  (2.12)  in  a  simple  manner,  but  there 
is  the  disadvantage  that  uncoupling  occ-ors  only  in  the  case  of  a  horizontal  mag- 
netic field  (i.e.,  when  y_  =  0  as  in  magnetic  equatorial  regions). 

An  alternate  set  of  equations  can  be  obtained  by  the  method  outlined 
at  the  end  of  Section  2,  The  motivation  for  the  particular  transformation  to  be 
used  is  implicit  in  the  interpretation  to  be  given.  TIlus  we  take  a  variable 
Lorentz  transformation  (i,e.,  rotation  through  a  complex  angle)  of  the  electric 
and  the  magnetic  field  components  separately.  However,  since  we  wish  the  new 
components  to  form  tv;o  orthogonal  pairs,  the  transformation  applied  to  the  mag- 
netic components  shall  be  the  inverse  of  that  applied  to  the  electric  components. 
Explicitly  we  introduce  the  new  variables  U(z)  as  in  equations  (2.8)  with 

U 

0 


(2.13)  U(z)  = 


^,(z) 
^2(0 
.-"2(2); 


.  P(z)= 


"JZJu) 


rr-1 


.P'MO= 


Jl^v.Hz) 


\  1  0  0 


Uj 


0 

1 


u 

1 

0 


0 

1 


\ 


u 


\  1  0  0 -y 


rr-l, 


The  quantity  u(z),  as  yet  unspecified,  is  determined  by  the  condition  that  P  AP 
Bhall  have  zero  off-diagonal  blocks  (i.e..  that  the  coefficient  matrix  in  the  equation 


-12- 


for  U(z)  "be  avich.   that  the  Tmcoupling  of  the  system  is  optimal  under  the  given  condit- 
ions), IVom  tthe  matrix  (2.10)  and  (2,13)  we  find  that  this  condition  reqiiires 


(2.14) 


u2  + 


[M^ 


+  1  =  0 


Tliis  equation  has  two  solutions  u  and  -•  »  where 


«  =  (^)  -  7(^)^-1 


(2.15) 


The  equation  satisfied  "by  Tl(z)  is  given  "by  equation  (2.9)  where,  with  the  above 
choice  of  u(z). 


0-100 


T      ■-^"(z)  0   0   0 
[2.16)  P"-^  AP  =  '   ^ 


0   0   0  xliz)!*"       ^^ 


u 


1-u 


0   0   10 


•   du    ♦ 
Here  we  have  introduced  the  notation  u  =  t~  and 

dz 


\?(z)     =     a+-    =     P-Va 
1  —  u    — 


\^(z)     =     a  +  Yu  =     P  -  Jf 


=     1  - 


-A 


0  0       0 

0  0       10 

0  10       0 

;-l  0 


0       0  J 


Tf 


T 


"-H  ^Tbyt/f^F^cby;  *^ 


This  system  of  equations  is  coupled  in  a  sli^tly  more  conrplicated  manner  than  the 
previous  system  "but   still  has  only  one  coupling  coefficient.    Eliminating  \^  ^  and 
iD  „  from  this  system  yields 


It  should  "be  pointed  out  that  +  \.  (z)  and  ±\^i7.)  are  the  eigenvalues  of  the  matrix 

A(z)  given  in  equation   (2.10).     The  sqtiares  of   these  eigenvalues  are  the  ordinary 
and  extraordinary  dielectric  constants,   e^  and  e_,  obtained  hy  Eydheck  jY)   and 

others. 


-1> 


(2.18) 


cl\(z) 


■'T,  ■  , 


2  o 

c 


te)1v^'=-&)s^-<J"= 


dz' 


C  1-U 


"a*') 


'-Ki^)3r-&(;^)''i 


These  are  Ryd"beck*s  [l]  well-knovm  equations  and  have  Tjeen  examined  by  variovis 
authors,  It  is  clear  that  the  coupling  enters  in  a  more  complicated  manner  here 
than  in  equations  (2.12).  However,  equations  (2,18)  uncoijple  not  only  in  the 
case  of  a  horizontal  magnetic  field,  hut  also  for  a  vertical  one,  since  then 
y_  =  0  and  hence  u  =  -i.  If  we  start  with  equations  (2,18)  then  we  mast  use  a 
limiting  process  for  u->  -i. 

Figure  1  represents  the  transformation  given  hy  equations  (2,8)  and  (2.1'?). 


^-^z 


^  >^l 


— >^ 
ii 

^2 


V, 


->•  X 


Figure  1 
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The  orthogonal  pair  of  field  components  V. ,  I  is  essentially  replaced  "by  an 
orthogonal  pair  tt  .  ^ .  whose  directions  are  related  co-unterclockwise  through 
an  angle  |6  from  the  directions  of  the  original  components,  A  similar  relation- 
ship holds  "between  the  remaining  pairs.  However,  it  should  he  recalled  that  the 
new  components  are  not  obtained  by  a  simple  rotation  of  the  vectors  involved. 


4,   Oblique  Incidence  (Vertical  magnetic  field) 

If  the  wave  is  normal  in  the  xz  -  plane  we  have  q.  =  0,  and  a  vertical 
magnetic  field  requires  y^  ~  7o  ~   Q«  ^®  coefficients  matrix  of  the  sjcstem  (2,6) 
becomes 

A 


(2.19) 


Mz) 


u  0  0 

0  0  V 

0  0  -p 

0  -1  0> 


where 


(2.20) 


-  P 


.     /<  =  1-  -^ 


(fi^-T^) 


Again  we  see  that  if  Y  were  zero  there  would  be  no  coupling.  If  we  eliminate 
I.  and  I  from  the  system, (2.6)  with  the  above  coefficient  raatrii  we  have 


(2.21) 


dV(z)    /;dV(z)      2 

=*= -  4 7 +  =::  >ua 

2     ^   dz 


dz 


dz^ 


/.a7^(z)  =  -^  A^Y^iz) 


+     ^       3  YM)     =    2-  V   V^(z)   . 
c  c 


Ihe  remaining  field  components  are 


I^(z)  =  p  ^^(z)   . 


V'>  '  'I 


p6     d7i(z) 
[(l.p2)^x]  "^^ 
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The  coupling  in  equations  (2,2l)  is  more  complicated  tiian  that  of  the  system  in  the 
previous  Section  "because  of  the  fsictoryu,  which  "becomes  unity  for  normal  incidence. 
Equations  (2,21)  are  equivalent  to  those  of  Wilkes  [5J  and  others.  This  may  "be 
seen  "by  multiplying  the  first  of  equations  (2.21)  "by  7^  aad  the  second  "by  (+  i) 
and  adding. 

The  method  of  Section  2  will  be  used  to  o"btain  an  equivalent  form  of 
equations  (2.2l),  A  comparison  of  (2,10)  with  (2,19)  suggests  a  transformation 
similar  to  (2,13).  However,  no  choice  of  u  greatly  simplifies  P  ~  AP.  Thus 
we  modify  the  transformation  "by  taking  U(z)  as  in  (2,13)  aad 


(2.22)  P(Z)  = 


h. 


./.V 


;  P''^(z)  = 


/vV 


r,-l 


We  pick  V(z)  such  that  P"  AP  is  simplified  as  much  as  pos8i"ble;  this  requires 


^*^' 


(see  [2],  Section  lO)  that 

(2.23)         V^*  (-^7^) 

This  equation  has  the  solution  T  and  rj  "7,  where 

\2Vy-/\2/iV/         /U 


=     1 


[sd-p^)-:^ 


^. 


[6(l-p^)-3 


6  -  X 
6(l-p^)-X 


> 


Now  the  coefficient  matrix  of  equation  (2,9)  has  the  components 


(2.25) 


and 


P"^  AP 


-1 

0 
0 
0 
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(2.25)  cont'd 


p-^  = 


lyC<v2 


V 


0  0  0  1 
0  0  -/>^  0 
0-100 
/^     0       0       0 


^ 


.  l-/^n 


Here  we  have  introduced 


(2.26) 


X^(z)  =  P  -  V/Y  =   //<  (a  +  W) 


X|(z)  =  P-V  T=   />(a+^) 


Note  that  again  +^.  (z)  and  i\Az)   are  the  eigenvalues  of  the  A(z)  now  given  by  (2.19). 
There  are  now  three  distinct  co-upling  factors: 


7 


1  -/AV' 


>U7 

1  -z^fv" 


1  -/x7' 


for  the  system  (2.9).  They  reduce  to  one  if  ^/X  =  1  and  then  we  have  the  case  of 
Section  3.   If  in  eq,  (2.9)  we  use  the  values  given  in  (2.25)  and  the  first  of 
(2.13)  and  then  eliminate  ^  .  and  ^  ,  we  ohtain 


dV(z) 


-[x|*a2-;-^h2*^]n, 

=    -     ac  +  c  -/Uab  -  ^^  I  n^  -   I c  +/^  "b    ri 


cL  TT,(z)  /U    cLtt  (z)        , 


Here  we  have  introduced 
(2.28)  a  = 


/^  r 


2(1  -yUV^)* 


h  = 


(1  -/x-T^)'     "  (1  -/Xt2) 


The  equations  (2.27),  while  equivalent  to  (2,2l),  are  by  no  means  simpler  and  seem 
to  have  only  the  virtue  that  they  reduce  to  the  I^dbeck  equations  (2.18)  when 
A*  =  1. 
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